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METHOD FOR CHARACTERIZING TUNABLE LASERS 



Related Applications 

This application claims priority from U.S. Provisional patent application 
60/41 1 ,858, filed on Sept 18, 2002, and which is Incorporated by reference. 

Field of the Invention 

The present invention is directed generally to the characterization of 
semiconductor lasers, and more particularly to approaches to characterize 
semiconductor lasers that are tunable. 

Background 

Tunable diode lasers have become widely accepted as important 
features of optical communications systems. Tunable lasers both simplify the 
maintenance of a dense wavelength division multiplexed (DWDM) 
communications system as well as enabling new network concepts. Some of 
the most important applications for tunable lasers include inventory control, 
frequency conversion, heterodyne detection, dynamic capacity allocation, and 
optical packet switching. 

Inventory Control 

The present light sources used in WDM optical transmitters are typically 
distributed feedback laser that emit light at a fixed frequency. As the number of 
channels in DWDM communications systems increases, for example to 80 or 
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more, carriers and system manufacturers are faced with the increased costs of 
maintaining a large inventory of spare transmitter laser. The availability of 
tunable lasers may greatly reduce the complexity of a DWDM transmitter unit 
by providing a single laser that may be programmed to emit over several, if not 
5 all, of the optical channels. A tunable laser may also be programmed to emit 
light at any desired optical frequency between the set optical channel 
frequencies. 

Frequency Conversion 

Some advanced optical communications system architectures require 

10 tunable frequency conversion. One approach to realize frequency conversion 
is to use a tunable laser as the transmitter in a transponder arrangement. In a 
transponder, the incoming optical signal is converted to an electrical signal by a 
photodiode, the electrical signal is amplified, as well as perhaps being 
reshaped to retimed, and is subsequently applied to an external modulator that 

15 modulates the output of the tunable laser. In future applications, all optical 
frequency conversion may be preferred, which will also require a tunable laser 
to generate the new carrier wave on which the data signal from the original 
carrier will be imposed. 

Heterodyne detection 

20 As the channel separation in DWDM systems decreases, it becomes 

harder to separate adjacent channels using optical filtering. With heterodyne 
detection, the problem of separating the densely packed frequency channels is 
moved to the electrical domain, where filters that are more selective are 
generally available. One element in a heterodyne receiver is a tunable laser, 

25 which acts as a local oscillator. The light from the local oscillator is combined 
with the light carrying the data, and the composite signal is detected using a 
pair of photodetectors. Thus, a radio frequency signal is generated having a 
frequency equal to the difference in frequency between the local oscillator and 
the data signal. To detect the signal at a particular frequency, the local 
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oscillator is tuned to within a few GHz of ttie carrier frequency of the signal. 
The local oscillator is then operated as a voltage controlled oscillator to keep 
the difference frequency constant. The neighboring channels may be 
suppressed by passing the composite signal through a band-pass filter 
5 centered at the difference frequency. 
Dynamic capacity allocation 

Due to the large growth of the Internet, with its rapidly shifting traffic 
patterns, networl< operators are seeking systems that allow rapid 
reconfiguration and dynamic capacity allocation. One approach to this problem 
10 includes a network with fixed frequency routers at its nodes, in which tunable 
lasers are used to dynamically set up paths across the network - changing 
frequency changes the physical path. 

Optical packet switching 

Although telecommunications traffic will soon be mostly dominated by 
1 5 data traffic, current optical network architectures do not take the "bursty" nature 
of this data traffic. Connections are typically set up between points for hours, 
days or longer. Data traffic is, however, transmitted in short packets, which are 
routed across the network without setting up long term point-to-point 
connections. Using rapidly tunable lasers, one can develop an optical packet 
20 network in which data packets are routed based on the frequency of their 
carrier waves. Such an architecture requires rapid tunability of the laser. 

In view of this wide range of applications, a broad spectrum of 
technological solutions for tunable lasers has been proposed: external cavity 
lasers; micro-electromechanically tuned, vertical cavity surface emitting lasers 
25 (MEMS-VCSEL); cascaded temperature-tuned distributed feedback (DFB) 
lasers; and various types of lasers based on the use of a distributed Bragg 
reflector (DBR). 

DBR-type lasers are particularly useful since the technology is the most 
mature, and the components are monolithic, which reduces the costs of 
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packaging and assembly, and permits the easy integration with additional 
components such as optical amplifiers or modulators. Moreover, DBR-type 
lasers can achieve the fast tuning times required for optical packet switching. 
There remains the problem, however, of characterizing a tunable laser, 
5 such as a DBR-type laser, after it has been fabricated. Correct characterization 
is required to ensure that the user knows how the laser wavelength tunes when 
one or more tuning parameters are changed, 



Summary of the Invention 

10 Generally, the present invention relates to an approach to characterizing 

tunable lasers. One particular embodiment of the Invention is directed to a 
method of characterizing a semiconductor laser having at least first and second 
tuning sections controlled by respective first and second tuning cunrents. The 
method includes measuring power output from the laser as a function of the 

15 first and second tuning currents, and creating an image of power as function of 
the two tuning currents. The image is analyzed to determine different modes, 
each mode corresponding to limited ranges of the first and second tuning 
currents. A preferred combination of the first and second tuning currents is 
determined for each mode and an acceptable operating region is defined for 

20 each mode. 

The above summary of the present invention is not intended to describe 
each illustrated embodiment or every implementation of the present invention. 
The figures and the detailed description which follow more particularly 
exemplify these embodiments. 

25 Brief Description of the Drawings 

The invention may be more completely understood in consideration of 
the following detailed description of various embodiments of the invention in 
connection with the accompanying drawings, in which: 
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FIG. 1 schematically illustrates a three section, distributed Bragg 
reflector laser; 

FIG. 2A schematically illustrates a sampled grating (SG) laser; 

FIG. 2B shows a graph of the reflectivity spectra of the two sampled 
5 gratings of the laser of FIG, 2A; 

FIG 3A schematically illustrates a grating coupler with rear sampled 
reflector (GCSR) laser; 

FIG. 38 shows a graph of the reflectivity spectrum of the sampled 
reflector and the transmission spectrum of the grating coupler of the laser 
10 illustrated in FIG. 3A; 

FIG. 4 schematically illustrates a tunable laser system including control 
electronics; 

FIG. 5 schematically illustrates a tunable laser characterization system; 

FIG. 6 schematically illustrates a laser system being tuned to particular 
15 frequencies based on a look-up table of laser characteristics generated using 
the characterization system illustrated in FIG. 5; 

FIG. 7 shows the side mode suppression ratio measured as a function of 
front and rear reflector currents for an SSG laser; 

FIG. 8 shows active selection voltage measured as a function of the front 
20 and rear reflector currents for an SSG laser; 

FIG. 9 shows a graph of the output from a frequency discriminating filter 
measured as a function of the coupler and reflector currents for a GCSR laser; 

FIG. 10 shows a graph of output power measured as a function of the 
coupler and reflector currents of a GCSR laser; 
25 FIG. 1 1 shows a graph of fiber-coupled output power and estimated 

frequency as a function of reflector current 

FIG. 12 shows a color-scale map of output power as a function of 
coupler current and reflector current, measured for decreasing reflector 
currents; 
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FIG. 13 shows a color-scale map of estimated frequency as a function of 
coupler and reflector current, measured for decreasing reflector currents; 

FIG. 14 shows a color-scale map of output power hysteresis as a 
function of coupler and reflector currents; 
5 FIG. 15 shows a pre-processed, color-scale map of output power as a 

function of coupler and reflector cun-ents; 

FIG. 16 shows a color-scale, segmented power image as a function of 
coupler and reflector currents; 

FIG. 17 shows a color-scale mode Image as a function of coupler and 
10 reflector currents, illustrating bands and columns; 

FIG. 18 shows a color-scale mode image, with ellipses fitted to the 
modes, as a function of coupler and reflector currents; 

FIG. 19 shows a graph of fiber-coupled output power and frequency as a 
function of phase current; 
15 FIG. 20 shows a graph of light power as a function of gain current for a 

single operating point; and 

FIG. 21 schematically illustrates determination of a mode boundary 
using a watershed technique. 

While the invention is amenable to various modifications and alternative 
20 forms, specifics thereof have been shown by way of example in the drawings 
and will be described in detail. It should be understood, however, that the 
intention is not to limit the invention to the particular embodiments described. 
On the contrary, the intention is to cover all modifications, equivalents, and 
alternatives falling within the spirit and scope of the invention as defined by the 
25 appended claims. 

Detailed Description 

In general, the present invention is directed to an approach for 
characterizing the operational characteristics of a tunable semiconductor diode 
laser based on the various tuning parameters. In particular, this approach 
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permits the user to define and characterize different areas of stable operation. 
The characterization approach may be used at different stages throughout the 
life of a laser. For example, the approach may be used both before and after 
initial burn-in, in order to determine the change in the laser's characteristics as 
5 a result of the burn-in process. 

As indicated above, lasers that use a distributed Bragg reflector (DBR) 
are particularly useful as tunable lasers for optical communications. DBR-type 
lasers typically use at least one, and often use two or more, different sections 
for tuning emission frequency of the laser. These different tuning sections are 

10 typically operated by injecting independently adjustable currents through the 
respective tuning sections. Therefore, two or more independently adjustable 
currents are often injected into a DBR-type laser for operation, a gain current 
injected into the active region of the laser to produce optical gain and one or 
more tuning currents to control the frequency of the light output form the laser. 

15 A schematic cross-section through one particular embodiment of a DBR 

laser 100 is presented in FIG. 1 . The laser includes three sections, the active 
section 102, the phase section, 104 and the reflector section 106. The 
reflector section 106 Includes a diffraction grating 108 having a period A that 
reflects light in a narrow band centered on the Bragg frequency, vb, given by: 

20 v»= 

' 2«(v,)A 

where n is the effective refractive index of the waveguide 1 10 in the 
reflector section 106. The waveguide 110 reaches between the reflector 
section 106 and the output coupler 1 12 at the other end of the active section 
102. Typically, the output coupler 1 12 is a cleaved facet of the semiconductor 
25 laser 100. The major laser output 1 14 is directed through the output coupler 
112. Electrodes 1 16, 1 18 and 120 are respectively disposed over the active 
section 102, the phase section 104 and the reflector section 106 for injecting 
independent currents into these sections 102, 104 and 106. A common 
electrode 122 is typically positioned under the laser substrate 124. 
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The effective refractive index, n, is changed by injecting current into the 
reflector section 106, resulting in a concomitant change in the Bragg frequency. 
The cavity mode is tuned to the Bragg frequency by adjusting the current 
passing through the phase section 104. The quasi-continuous tuning range of 
5 the DBR laser 100 is limited by the maximum change in the effective refractive 
index of the waveguide 110 that can be accomplished by current injection. 
Through careful optimization of the waveguide structure, tuning ranges up to 2 
THz may be achieved. This is, however, significantly less than the gain 
bandwidth of the active section 102, which may reach about 10 THz for InP 

1 0 semiconductor lasers, a type of semiconductor laser often used to generate 
light used for optical communications, for example in the wavelength band 1500 
nm - 1 620 nm. This tuning range is also less than the gain bandwidth of the 
erbium-doped fiber amplifier (EDFA), an amplifier that finds widespread use in 
long-haul optical communications systems. 

15 An example of another type of laser based on a DBR is a sampled 

grating (SG) laser 200, an embodiment of which is schematically illustrated in 
FIG. 2A. This laser 200 includes first and second reflector sections 202 and 
204, with an active section 206 and a phase section 208 disposed between the 
reflector sections 202 and 204. A waveguide 210 guides the light through the 

20 active and phase sections 206 and 208, between the reflector sections 202 and 
204. Electrodes 212, 214, 216 and 218 are disposed above respective 
sections 202, 204, 206 and 208 to permit the injection of current into the 
different sections independently. A common electrode 220 is typically 
positioned under the laser substrate 222. 

25 In the illustrated embodiment, each reflector section 202 and 204 

includes a sampled diffraction grating, each of which has a comb-shaped 
reflection spectrum. The reflectivity spectrum of a grating can, within certain 
limits, be described by a Fourier transform of the grating function. As is known 
from Fourier theory, a periodic modulation of a carrier wave yields a comb- 
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shaped spectrum centered on the frequency of the carrier wave, having a peal< 
separation equal to the frequency of the modulation. Therefore, a grating that 
is modulated with a period As manifests a reflectivity having reflection peaks at 
a set of frequencies, Vk, given by: 



2nK) 



where k Is an integer value. 

The periodic modulation may be in amplitude, as in a sampled grating, or 

10 in phase. In the latter case, the structure us commonly called a super-structure 
grating (SSG). One advantage of the SSG is that there is more freedom in 
designing the relative strengths of the various reflectivity peaks, although this 
comes at the cost of more complex semiconductor processing. 

An SG-DBR laser, or SSG-DBR laser is tuned using a Vernier effect, as 

15 is illustrated with reference to FIG. 2B. The reflector sections 202 and 204 are 
designed such that the reflectivity peak separation (6f) of the front reflector 
section 202 and the reflectivity peak separation (6r) of the rear reflector section 
204 are slightly different: 5f - 6r = A5. This is achieved by providing the 
sampled sections of grating 224 and 226 with slightly different modulation 

20 periods, Asf and Asr respectively. The laser operates at that frequency of the 
cavity mode that falls within the frequency band where a reflectivity peak of the 
front reflector section 202 overlaps with a reflectivity peak of the rear reflector 
section 204. The cun-ent through the phase section 208 may be adjusted to 
precisely align the cavity mode with coinciding reflectivity peaks, and to tune the 

25 laser to the desired optical channel frequency. If one of the front or rear 
reflector sections 202 or 204 is tuned by A6 or more, then two neighboring 
reflectivity peaks coincide and, with proper adjustment of the current through 
the phase section 208, the frequency of the light emitted changes by 6f » A6 
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where the front reflector current is changed and by - 6r » A5 where the rear 
reflector current is changed. 

It will be appreciated that, instead of changing the tuning current of only 
one reflector section 202 or 204, the laser 200 may also be tuned by adjusting 
the tuning current of both reflector sections 202 and 204. The SG-DBR or 
SSG-DBR laser is described in greater detail in U.S. Patent No. 4,896,325, 
incorporated herein by reference. 

An embodiment of a grating-assisted coupler, sampled reflector (GCSR) 
laser 300, is presented in FIG. 3A. The laser 300 includes four sections, a gain 
section 302, a coupler section 304, a phase section 306 and a reflector section 
308, each typically integrated on the same substrate 310. The gain section 302 
includes an active waveguide 312, and may include a quantum well structure to 
provide optical gain. A gain electrode 31 1 is disposed over the gain section 
302 to permit injection of cun'ent through the gain section 302. A common 
electrode 313 is typically disposed over the bottom surface of the substrate 
310. 

The waveguide 312 extends into the coupler section 304 as a first 
waveguide 316. A second waveguide 318 lies close to the first waveguide 316. 
A grating structure 320 is disposed near the second waveguide 318. The 
grating structure 320 is illustrated above the second waveguide 318, but may 
optionally also be placed between the first and second waveguides 316 and 
318, or below the first waveguide 316. A coupler electrode 322 may be 
disposed over the coupler section 304 to permit injection of current through the 
coupler section 304. 

The second waveguide 318 couples to a phase waveguide 323 through 
the phase section 306 and into the reflector section 308. A phase electrode 
324 may be disposed over the phase section 306 to permit injection of cun-ent 
through the phase section 306. 
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The reflector section 308 includes a reflector structure 326 disposed 
near the reflector waveguide 325 that is coupled to receive light from the phase 
waveguide 323. In the illustrated embodiment, the reflector structure is a 
sampled Bragg reflector, although the reflector structure 326 may be any type 
5 of reflector structure that provides the desired reflective characteristics. A 

reflector electrode 328 may be disposed over the reflector section 308 to permit 
injection of current through the reflector section 308. 

The GCSR laser 300 is able to produce light in a single longitudinal 
mode that is widely tunable over a large wavelength range, and is particularly 

1 0 suitable for use as a source in dense wavelength multiplexed (DWDM) optical 
communications systems. The laser cavity is formed between the output facet 
330 and the reflector section 308. In other embodiments, not Illustrated, the 
output coupler of the laser 300 may be a wideband Bragg reflector, as is 
described further in U.S. Patent Application Serial No. 09/915,046, incorporated 

1 5 herein by reference. The use of a grating-assisted coupler and sampled Bragg 
reflector for tuning a laser is described further in U.S. Patent No. 5,621 ,828, 
incorporated herein by reference. 

By using a grating structure, having a period Ac, efficient power transfer 
may be effected between the waveguides 316 and 318 over a limited frequency 

20 band around the coupling frequency vc, given by: 

c 

^•'~AjnR(vJ-ns(vJ] 

where nR (vc) and ns (vc) are the effective refractive indices for light at 

frequency Vc for the waveguides 316 and 318 respectively. Since the coupling 

frequency depends on an index difference, a small change in either nR or ns 

25 yields a frequency change given by : 

Av^^AK-nJ 
Vc HR-ns 
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The coupler section 304 by itself may not provide sufficient frequency 
selection to achieve single mode operation with good side mode suppression, 
and so a sampled reflector section 308 may be used to provide a reflectivity 
spectrum that includes a number of highly reflection peaks 356, illustrated in 
5 FIG. 3B, separated by regions of wavelength where the reflectivity Is low. In the 
particular embodiment illustrated in FIG. 3B, the separation between the 
different reflection peaks 356 is AAp. The coupler section 304 has a relatively 
broad transmission spectrum 358, which is wavelength tunable by injecting 
different amounts of current via the coupler electrode 322. Therefore, the 

1 0 transmission window 358 of the coupler section 304 may be tuned to select a 
single reflection peak 356 of the reflector section 308, thus selecting a single 
longitudinal mode for oscillation. Since the reflectivity peaks 356 of the reflector 
section 308 are also wavelength tunable by injecting different amounts of 
cun-ent through the reflector electrode 328, the laser 300 may be made to 

15 oscillate on a single longitudinal mode at substantially any selected wavelength 
within the operating wavelength range. The oscillating wavelength may be fine- 
tuned by adjusting the current injected through the phase section 306 via the 
phase electrode 324. 

For a laser used as a source In optical communications having a 

20 wavelength in the range 1 500 - 1 620 nm, a typical wavelength range for long- 
haul fiberoptic communications, the lasers 100, 200 and 300 may be based on 
indium phosphide (InP), having an InP substrate. The waveguides 110, 210, 
316, 318, 323 and 325 are typically formed of a material having a higher 
refractive index than the surrounding material, in order to provide optical 

25 confinement. The waveguides 31 6, 31 8 323 and 325 may be, for example, 
formed from an Indium gallium arsenide phosphide (InGaAsP) alloy. The 
grating structure 320 may also be fomied from islands 330 of high refractive 
material, for example InGaAsP, spaced apart in a repetitive pattern. 
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Tuning a laser, that has different tuning sections controlled by different 
tuning currents, to a particular frequency with high side-mode suppression ratio 
(SMSR) may require simultaneous adjustment of up to three or more different 
tuning currents. Typically, the requirement for SMSR is that the adjacent mode 
5 be suppressed by more than 35 dB. 

Due to fabrication tolerances, the set of currents selected to produce 
light for a particular frequency channel may vary from laser to laser. For 
practical applications, the laser is, therefore, supplied with control electronics 
that contain a channel look-up table, for example stored in an EEPROM. An 

1 0 embodiment of a tunable laser system 400 is illustrated a block schematic 
diagram presented in FIG. 4. Such a laser system 400 may be incorporated in 
a DWDM transmitter unit. The laser 402 generates an output light beam 404, a 
portion of which may be directed to a wavelength detector unit 406, which 
generates an output signal 408 determined by the wavelength of the light in the 

15 light beam 404. 

A residual output beam 410, passing from the wavelength detector unit 
406, may carry optical output power not used in the determination of the 
wavelength. The residual output beam 410 may be used as the useful optical 
output from the laser 402. Where the output light beam 404 carries the main 

20 optical output from the laser 402, the wavelength detector unit 406 

advantageously uses only a small fraction, for example a few percent, of the 
output light beam 404, in order to increase the power in the residual output 
beam 410. 

A wavelength analyzer unit 412 may receive and analyze the output 
25 signal 408 from the wavelength detector unit 406 to determine the wavelength 
of the light beam 404. The analyzer 41 2 typically generates an error signal 414 
that is directed to a wavelength controller. The size of the error signal typically 
indicates the amount by which the measured wavelength of the laser deviates 
from a desired value. The error signal 414 is directed to a wavelength tuning 
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controller 41 6 that is connected to the laser 402 and controls the operating 
wavelength of the laser 402. The wavelength tuning controller 416 may, for 
example, direct different tuning currents to different sections of the laser 402. 
The wavelength tuning controller 416 may be incorporated with a laser 
5 controller 41 8 that includes the power supply 420 for providing power to the 
laser 402 and a temperature controller 422 that controls the temperature of the 
laser 402. The laser 402 may be coupled, for example, to a thermoelectric 
device 424 or other type of device for adjusting temperature. 

The wavelength tuning controller 416 may include a memory device 428, 

10 such as an EEPROM, that contains the lool<-up table that indicates the different 
tuning currents that are used to achieve a laser output at a particular channel 
frequency. The wavelength tuning controller 416 may also contain circuitry that 
provides compensation for the tuning currents applied to the laser, for example 
to compensate for drift in laser temperature, aging of the laser, or other effects 

15 that may change the optimum values of the tuning currents. Such current 
compensation may, for example, be based on the size of the error signal 
received from the wavelength analyzer unit 412. The wavelength tuning 
controller 416 may also be coupled to receive an external control signal 430 
that controls the optical channel on which the laser oscillates. The external 

20 control signal 430 may be received, for example, from an optical 

communications system controller that controls operation of the optical 
communications system. 

The laser 402 and wavelength detector unit 406 may be enclosed within 
a housing 426 to prevent environmental effects from affecting the operation of 

25 the laser 402 and the wavelength detector unit 406. The device 424 for 

adjusting operating temperature may also be located within the housing 426. 

In practice, the measurements required to characterize a laser should 
take as short a time as possible. Characterization of the laser includes 
measuring its tunability as a function of the different tuning currents and 
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showing that it can achieve desirable levels of SMSR. The measurements are 
preferably made within a time of a few minutes or less. The procedure is 
mostly concerned with determining the optical performance of the laser, for 
example, frequency tuning range, output power, side-mode suppression ratio, 
5 threshold and other performance related parameters. Several parameters 
(scalars) are extracted from measurement data and compared to the limits, 
typically given as maximum and/or minimum values, listed in the engineering 
specification for the laser chip. If the laser fails to meet these requirements at 
any time, the laser may be rejected. The procedure may be performed prior to 

10 burn-in and after each burn-in step, thus permitting the degradation in the 
performance of the laser to be monitored during the burn-in process. 

Steps of one particular embodiment of a process for characterizing a 
tunable laser are listed in FIG. 5. A characterization process need not include 
all of these steps, or may include additional steps. The first step 502 includes 

15 scanning the current of a first tuning element, typically not a phase tuning 
section, and measuring the resulting output power. The second step 504 
includes setting the laser to a tuning current that produces a relatively high level 
of output power. For example the tuning current may be set to the level 
associated with the maximum output power. The laser may then be aligned to 

20 a fiber for coupling to diagnostic equipment. A second scan of the tuning 
current is made while measuring output power and wavelength. 

The tuning current of a second tuning element may then be scanned, in 
step 506 to make some initial measurements of power and wavelength. Next, 
at step 508, the tuning currents of the first and second tuning elements are both 

25 scanned, and the power and wavelength mapped for the two dimensional 

tuning current space. The data obtained in step 508 are then analyzed at step 
510, to determine those combinations of tuning currents that result in stable 
single mode operation. Next, at step 512, the SMSR is measured for different 
wavelengths to obtain different operating points where the SMSR is reduced. 
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Then, at step 514, the frequency and output power is measured for the different 
operating points by tuning the phase section. Finally, the threshold is 
measured at step 516 for the different operating points. 

These steps need not all be performed in the same order in which they 

5 are listed, and some steps may be omitted. 

One embodiment of a system 600 that may be used for characterizing a 
laser is illustrated schematically in FIG. 6. The characterization of the laser is 
typically performed once the laser has been mounted on a carrier, or submount, 
and is described as being at the laser on carrier (LoC) level. The laser carrier 

10 602 is mounted on a probe mount 603 to make electrical contact with the 
various electrodes of the laser 604. A current controller 606 supplies drive 
current and tuning currents to the laser 604. For example, the current controller 
608 may supply a drive current, Id, and three tuning currents, 11, 12, and 13. 
More or fewer tuning currents may also be supplied. The light output from the 

15 laser 604 may be measured using a calibrated power photodiode 608, for 

example as may be used for step 502. The output from the laser 604 may also 
be directed to a power/wavelength measuring unit (PWU) 610, which may be 
fiber-based. The PWU 610 measures the power produced from the laser 604 
as a function of wavelength, as is discussed below. The PWU 610 is coupled 

20 to a digitizer/processor 612 that analyzes the data produced by the PWU 61 0. 
The digitizer/processor 612 also controls the operation of the current controller 
606 so that the data obtained from the PWU 610 may be related to the 
associated values of the tuning currents. The digitizer/processor 612 may be 
programmed to run the characterization program automatically. An output 

25 device 614, for example a printer, screen, and/or the like, coupled to the 
permits the user to view the results of the characterization process. 

Two different embodiments of PWU are illustrated in FIGs. 7A and 8A. 
In the first embodiment of PWU 700, illustrated in FIG. 7A, incoming light 702 
from the laser is split by a beamsplitter 704 into two beams 706 and 708. The 
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first beam 706 is directed to a first photodetector 712, sucli as a photodiode, to 
monitor power. Tlie second beam 708 is directed through a filter 714 having a 
known transmission characteristic to a second photodetector 716. The 
transmission of the filter 714 increases or decreases with wavelength. The 

5 ratio of the signals produced by the two photodiodes 71 2 and 71 6 permit an 
estimation of the wavelength or, equivalently, the frequency, of the light 702. 
The graph illustrated in FIG. 7B shows a characteristic plot of photodetector 
signals as a function of frequency. The signal 722 generated by the first 
photodiode 712 is independent of frequency while the signal 726 generated by 

1 0 the second photodiode 71 6 is frequency dependent. 

In the second embodiment of PWU 800, illustrated in FIG. 8A, incoming 
light 802 from the laser is split by a beamsplitter 804 into two beams 806 and 
808. The first beam 806 is directed to a first photodetector 812 to monitor 
power. The second beam 808 is directed to a filter 814 having a known 

15 transmission characteristic. The filter splits the beam 808 into two beams 816 
and 818. The beam 816 is directed to a second photodetector 820 to monitor 
the amount of light reflected by the filter 814. The beam 818 is directed to a 
third photodetector to monitor the amount of light transmitted through the filter 
814. The transmission of the filter 814 increases or decreases with 

20 wavelength, and so the amount of light in beam 816 increases as the amount of 
light in beam 818 decreases, and vice versa. The graph illustrated in FIG. 8B 
shows a characteristic plot of photodetector signals as a function of frequency. 
The signal 832 generated by the first photodetector 812 is independent of 
'frequency while the signals 840 and 842, derived from photodetectors 820 and 

25 822 respectively, show a dependence on the frequency. The PWU 800 

generally permits a more accurate estimation of the frequency of the Incoming 
light than the PWU 700. An advantage of the PWU 800 is that it allows a more 
accurate estimation of the optical frequency of light input to the device. 
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The following characterization procedure is described in terms of its 
application to a GCSR laser. It will be appreciated, however, that the 
procedure may also be applied to other types of tunable laser, such as the 
other types of tunable laser discussed above, including the SG-DBR laser and 
5 the DBR laser. In the case of the GCSR laser, the currents applied to the laser 
are: the gain current, Ig, applied to the gain section 302, the coupler current, Ic, 
applied to the coupler section 304, the phase current, Ip, applied to the phase 
section 306, and the reflector current, Ir, applied to the reflector section 308. 
According to one embodiment of the invention, step 502, the first scan of 

10 the first tuning element, includes a scan of the coupler current, Ic. Under this 
step, the gain current, Ig, phase current, Ip, and reflector current, Ir, are each set 
to respective fixed values. For example, Ig may be set to a maximum specified 
value, while Ip is set to zero and Ir is set to a minimum specified value. 

The coupler current, Ic, is then swept over a range of values, for example 

15 0-40 mA, The output power from the laser is measured as a function of 
coupler current. An example of the result of such a measurement is shown in 
FIG. 9, which shows a plot of output power P(lc) as a function of Ic. The results 
may be analyzed by calculating a running average of the P(lc) curve. Also, the 
maximum output power, (P_max) and the corresponding coupler current 

20 (Ic-P^max) may be calculated. In the example illustrated in FIG. 9, P_max = 6.5 
mW and lc.P_max = 8 mA. 

According to one embodiment of step 504, the first tuning element is 
then scanned a second time, and the laser may be set with the same current 
values as in step 502, with the coupler current U set to lc.P_max. The output 

25 from the laser may then be aligned through an optical fiber. Alignment of the 
optical fiber to the laser at the maximum output power level reduces the 
difficulty in making the alignment. Measurement of the maximum power 
through the fiber permits a calculation of the fiber-coupling efficiency. 
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The coupler current, U, may then be swept while the output from the 
laser Is measured for both optical power and frequency. An example of the 
results of such a measurement is illustrated in FIG. 10, which shows both the 
fiber-coupled power, curve 1002, and the estimated frequency, curve 1004, as 

5 a function of U- This permits the operator to determine the range of coupler 
currents required to cover the desired frequency range. For example, where 
the desired frequency range of the laser is 192 THz - 1 96 THz, then the curve 
in FIG. 10 shows that the range of currents required to achieve this range is 
approximately lc_min = 7 mA and Ic^max = 26 mA. to make sure that the 

10 selected Ic is suffciently broad to cover the desired tuning range, the minimum 
(maximum) value of Ic may be decreased (increased) by some fraction, such as 
30%. A coupler current operating point, URScan may then be selected, for 
example by selecting a point that lies approximately in the middle of one of the 
stairs in the staircase-lil<e frequency v. Ic curve, curve 1004, preferably close to 

1 5 the maximum output power. 

The next step, step 506, is to scan the second tuning element which, in 
this particular embodiment, is. the reflector section of the GCSR laser. With Ig 
and Ip still at the same value as before, the coupler current is set to IcRscan. 
The reflector current, Ir, is then scanned from the minimum value to a maximum 

20 value. In the particular example, Ir is scanned from 0 mA - 40 mA. The output 
power and the laser frequency are measured as a function of If. An example of 
the results of such a measurement are illustrated in FIG. 1 1 , which shows 
output power, curve 1 102, generally sloping from about 2.75 mW at Ir = 0 mA to 
about 1 .2 mW at Ir = 40 mA. The other curve is the estimate frequency, curve 

25 1 1 04, plotted as a function of Ir. 

A value of lr_max is determined as that value of power required to tune 
the laser to the same frequency as the minimum value of Ir. In the particular 
example illustrated in FIG. 1 1 . Ir_max is about 25 mA. To ensure that the 
selected cun-ent range is sufficiently broad, the value of lr_max calculated from 
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the measurement data may be increased by a selected margin, for example 
30%. 

The next step, step 508, includes scanning both the first and second 
tuning elements. The gain and phase currents may be held at the same values 

5 as before. One of the currents, for example Ic, may be swept from Ic^mln to 
lc_max in a given number of steps. In the illustrated example, the given 
number of steps is 400. For each value of Ic, the reflector current, Ir, may be 
swept from 0 mA to lr_max and then back from lr_max to 0 mA, in a given 
number of steps. In the Illustrated example, the number of steps is 750. The 

1 0 optical power and estimated laser frequency may be measured for each 
combination of Ic and Ir. Various images representing these measurements 
may be made. One example of an image is presented in FIG. 12, which shows 
output power (color-coded) as a function of Ir (x-axis) and Ic (y-axis), measured 
for decreasing reflector curents. The color red represents relatively high power 

15 and the color blue represents relatively low power. The data are presented 
such that the upper left corner corresponds to the minimum values of U and Ir. 
Another example of an image may be formed for measurements taken when 
the value of Ir is increasing. Another example of an image may be formed from 
the hysteresis of the output cun-ents, for example the difference between the 

20 values of output power for increasing and decreasing values of Ir. 

Another example an image is an estimate frequency as a function of U 
and Ir, measured for either decreasing values of Ir, or increasing values of Ir. 
FIG. 13 shows estimated frequency for decreasing values of Ir. Higher 
frequencies are shown as red and lower frequencies shown as blue. The large 

25 frequency changes that occur when tuning the cun-ents are easily recognized, 
and con-espond to the laser frequency hopping from one cavity mode to 
another, commonly referred to as mode-hopping. 

The different power values may be scaled by dividing by the fiber 
coupling efficiency measured at step 504. 
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The next step, step 510, is to analyze the data taken in step 508. One 
approach to this is to form a hysteresis image. The difference between the 
power values for increasing and decreasing values of Ir, may undergo a 
thresholding process to eliminate small uncertainties in the measurements. 
One example of a thresholding process is to assign a pixel value of 1 to all 
points where the output power measured for increasing values of Ir is different 
from that measured for decreasing values of Ir by a given amount, for example 
5%. A pixel value of zero may be given to those pixels where the difference in 
power is less than the given amount. An example of a thresholded hysteresis 
image is presented in FIG. 14, in which only differences of more than 5% 
between power for increasing and decreasing values of Ir are shown. In this 
image, red corresponds to a large difference, while blue corresponds to a small 
difference. 

A frequency gradient may be calculated as follows: first, calculate a 
pseudo-Gaussian convolution kernel: 



u(i) = 
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exp 
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where £ is a small number (e.g. 0.0001). 
Next, derivative of the Gaussian convolution kernel is calculated: 

/ -.2 \ 



v(i) = — rexp 



2a^ 



where -N<i<N. 
These two kernels may then be convoluted: 
w(i) = 2:u(kMi-k) 
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The frequency image is convoluted with this kernel w(i), both in the x and y- 
directions. The results are then squared, added and the square root taken. 

The frequency gradient image may then be normalized by dividing the 
gradient values with the frequency separation between two neighbouring 
5 reflectivity peaks of the reflector, which should be equal to the frequency 
difference between two adjacent "bands" in the frequency image of FIG. 13. 

All values below a certain threshold value, for example, 0.2, may be set 
to zero in order to remove noise. All values above 1 may be set to one. The 
power image in FIG. 12 may then be multiplied by (1- the normalized frequency 
10 gradient image). This effectively lowers the image intensity in areas with high 
frequency gradient. An example of the resultant image is presented in FIG. 15. 

The processed power image shown in FIG. 15 may then be further 
analyzed, for example using a modified watershed algorithm, for example as 
discussed in "Watershed segmentation of binary images using distance 
15 transofrmations" Orbert, Bengstsson and Nordin, Proceedings of SPIE 

conferecne on Image Processing: Nonlinear Image Processing IV, San Jose, 
California, 1993; and "Watersheds in digital spaces: an efficient algorithm 
based on immersion simulations", Vincent and Soille, IEEE Transactions on 
Pattern Analysis and Machine Intelligence, vol. 13, pp. 583-598, 1991, both of 
20 which are incorporated by reference. 

The watershed algorithm is described briefly with reference to FIG. 21 , 
which shows a stylized cross-section through a power plot, for example as 
illustrated in FIG. 12. The watershed algorithm finds the boundaries between 
different modes by examining the gradient of the power curve. For example, 
25 the algorithm examines the gradient of the power curve around the current Iq. 
Since the gradient of the curve is different on either side of lo, the current lo is 
determined to be at a mode boundary. 

Since the frequency of the laser changes upon passing through a mode 
boundary, it is possible to verify the presence of a mode boundary, as 



Attorney Docket: 980.1131USU1 



determined using the watershed algorithm, by ensuring that the frequency also 
changes at the mode boundary current. Problems may occur if there is noise 
on the power curve. For example, the algorithm may assume that a local noise 
minimum is a mode boundary. Verification of a mode boundary using the 
5 frequency data reduces the possibility that the algorithm mis-characterizes 
noise as a mode boundary. Another possibility is that a power peak for a 
particular mode is not very high, and is assumed by the algorithm to be noise. 
Again, verification by comparing with frequency data may help to reduce the 
possibility that the algorithm fails to recognize a mode boundary. 

10 After application of the watershed algorithm, the segments in the image 

are sorted with respect to the frequency of the geometric midpoint . Segments 
are a power value at the midpoint (see, for example FIG. 12) that are less than 
some fraction of the maximum power, for example, 20%, may be removed. 
The result is presented in FIG. 16. 

15 The different, isolated segments in FIG. 16 represent different 

longitudinal modes of the laser. These modes may be sorted. First, the 
segments that touch the edges of the area may be removed. These segments 
represent modes that cannot be completely accessed using just Ic and Ir alone. 
For further processing, the segmented image may be divided into a 

20 number of vertical fields, for example 5 fields. The horizontal axis, the Ir axis, 
may be divided into parts with increasing width. The width of the different parts 
may increase linearly across the current range. This is based on the 
observation that the segments increase in width along the horizontal axis. A 
segment is said to belong to a certain field when its (geometric) midpoint lies 

25 between the left and right boundary of that field. 

For each field, the maximum area of a segment is first determined. 
Then, the average area is calculated of all segments that have an area 
between 20% and 90% of the maximum area. In other words, extremes are 
disregarded. Subsequently, those segments that have an area that is less than 
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25% of this average are removed. In this way, the small segments that lie 
squeezed between the larger segments in FIG. 10 may be removed. The 
remaining segments may then be sorted into bands based on the minimum and 
maximum y-coordinates and the y-coordinate of the midpoint. 

At the end, the bands of the different fields are connected to each other. 
The segments that remain after this process correspond to areas in which the 
laser operates in a single cavity-mode. These segments are, therefore, referred 
to as modes. 

The modes may also be divided into "columns", that is continuous lines 
may be laid over the modes, to connect vertically adjacent modes. The lines 
are shown as dotted lines 1702 in Fig. 17. This eases detection of any modes 
that may be missing from one of the bands. The resulting image is shown in 
FIG. 17. Bands are shown connected by dashed lines. 

A "workspace", generally an elliptic area, may be calculated for each 
mode. The workspace corresponds to a well-defined operating region that fits 
within the boundaries of each mode. To find the workspace of a mode, first, 
the mid-point of each mode (Xc, yc) is calculated, where: 



The sums are over all pixels of the mode, where n is the number of pixels. 
The elements of moment of inertia tensor, for axes through the midpoints may 
then be calculated: 





N 

ixx=E(yk-yc) 



k=1 



N 

lyy=Z(X|c-Xc)' 




From this, the principal moments of inertia may be calculated as: 




The slopes of the principal axes of the mode may be calculated as: 
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For each mode, the largest ellipse that satisfies the following conditions 
is calculated: 

a) the midpoint is (Xc.yJ; 

b) the ratio of the length of the principal axes of the ellipse is: 



c) the principal axis with length a lies along the line with slope m, ; and 

d) the entire ellipse lies within the mode. 

The ellipses for the modes in FIG. 17 are shown in FIG. 18. The size of the 
ellipses may be used as a criterion for selecting whether a laser is useful or not. 
For example, a laser may be rejected where more than a certain number of 
modes have ellipses with minor axes that are less than a particular threshold 
value. Such a characteristic may indicate that it will be difficult to obtain stable 
operation of such a laser, and the laser may be rejected. 

Typical output data for each mode are listed in Table I. 

Table I Output Data for Each Mode 



1 . Midpoint of the mode 

2. Band index. 

3. Column index. 

4. Gain current for the operation point (the gain current at 
which the image data were measured, i.e. the maximum 
specified gain current). 

5. Coupler current for the operation point (corresponding to 
the coordinate yj. 

6. Reflector current for the operation point (corresponding to 
the coordinate ). 
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7. Phase current for the operation point (the phase current at 
which the image data was measured, in this case 0 mA). 

8. Output power at the operation point (for example from FIG. 
12). 

5 9. Estimated frequency at the operation point (for example 

from FIG. 13). 

10. Area of the mode (in mA^). 

1 1 . Fraction of the area of the mode that shows hysteresis 
(calculated by overlaying FIG. 17 with FIG. 14). 

10 12. Relative size of the workspace in the reflector direction, 

defined as the ratio of the width of the ellipse along a 
horizontal line through the center to the mode separation in 
the horizontal direction. See the definition of Axc below. 

1 3. Relative size of the worl<space in the coupler direction, 
1 5 defined as the ratio of the height of the ellipse along a 

vertical line through the center to the mode separation in 
the vertical direction. See the definition of Aye below. 

14. Ellipse parameters (size of the principal axes and slope). 

20 The mode separation Ax^ in the horizontal direction may be calculated 

as: 

If both the previous (/c-1 ) and the next (/c+1 ) 
mode in the band exist. 
If only the next mode in the band (/c+1) exists. 
If only the previous mode in the band (/c-1) 
exists. 



^c,k ~ ^c,k ^c,k-1 



The mode separation Ay^ in the vertical direction may be calculated as: 

1 / \ If both the previous (/-I ) and the next (/+1 ) 

Ay.,=^(y.,.-y.H) column exist. 

If only the next mode in the column (/+1) 

^Ycj =yc.H -Ycj 

exists. 

If only the previous mode in the column (/-I ) 

^Vcj yc.j ^ yc,j-i 

exists. 
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Possible errors in the image analysis may be reported in a separate error list. 

The next step, step 512, is to measure the frequency and the side mode 
suppression ratio (SIVISR) for the different modes. For each of the operating 
points determined in step 510, the following parameters are measured: 
5 frequency (nu, usually measured in THz); 

side mode suppression ratio (in dB) 

frequency of the strongest side mode (nu_sm) 

These data are then analyzed, for example in the following manner. 
First, the average mode separation of a cavity that includes only the gain, 
10 coupler and phase sections of the GCSR laser, also referred to as the GCP 
mode separation, is calculated. This is given by the average difference in 
frequency between two operation points that lie in the same band, in other 
words have the same band index, and are adjacent to each other, in other 
words, whose column indices differ by 1. This average mode separation may 
15 be referred to as Deltanu_GCPMode . 

The average peak reflector separation, Deltanu_ReflPeak, is then 
calculated as the average difference in frequency between two operation points 
that lie in the same column and are adjacent to each other, in other words have 
a band index that differs by 1 . 
20 The coupler current, lc.nu_min, corresponding to the lowest frequency of 

the required tuning band, nu_min, is then calculated. This is done by finding 
two neighbouring operation points whose frequencies straddle nu_min. The 
value of lc.nu_min is calculated by linearly interpolating between the Ic values 
for the selected operation points. 
25 Next, the coupler current, lc.nu_max, that corresponds to the highest 

frequency of the required tuning band, nu_max, is calculated. This is done by 
finding two neighbouring operation points whose frequencies straddle nu_max. 
The value of lc.nu_max by then be calculated by linearly interpolating between 
the Ic values for the operation points. 
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The coupler current, lc_max, that corresponds to nujimit = nu_max + 
0.75*Deltanu_ReflPeak + Deltanu_GCPMode may then be calculated. This 
corresponds to the maximum coupler current needed to be able to measure a 
mode plane image that contains all modes (in their entirety) needed to cover 
5 the desired frequency tuning range from nu_min to nu_max. This is done by 
finding two neighboring operation points having frequencies that straddle 
nujimit. The value of lc_max may then be calculated by interpolating linearly 
between the Ic values for the selected operation points. 

If no two operation points that straddle nujimit can be found, then a 

10 linear extrapolation technique may be used to calculate the value of lc_max. 
For example, a straight line may be fitted to the curve of obtained when plotting 
frequency as a function of coupler current, and a value for lc_max may be 
extrapolated by extending the line to nujimit. 

The average reflector current, lr_max, needed to tune the reflector by 

15 the peak separation, starting from the minimum specified reflector current may 
then be calculated. This may be done by finding the first operation point 
(OP[band][col]) that has Ir > lr_min, for each band except the first band. From 
this OP and the previous OP (OP[band][col-1]), the start frequency for the band 
nu_start may be calculated by interpolating linearly between the frequencies of 

20 the two operation points. In other words, the frequency at lr_min is calculated 
as if the frequency increases linearly between the Ir-values of the two OP. If 
there is no previous OP, then the next OP (OP[band][col+1]) is taken and an 
extrapolation back to lr_min is made. If neither the previous nor the next OP 
exists, the band may be ignored. 

25 Subsequently, the two neighbouring OP are found within the band that 

have frequencies that straddle than nu_start + Deltanu^ReflPeak + 
Deltanu^GCPMode. The value of lr_max may be calculated for this band by 
interpolating linearly between the Ir values for the neighboring operation points. 
The average average lr_max may then be calculated across all bands. 
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The coupler tuning efficiency, TuningEfficiencyCoupler, typically 
measured in THz/mA, may be calculated at lc.nu_min. This is done by 
considering the curve obtained when plotting frequency as a function of coupler 
current, Ic, for all operation points. A straight line is fitted to this curve at 
5 Ic = lc.nu_min, using any suitable line fitting technique. The value of 
TuningEfficiencyCoupler is the slope of the fitted line. 

The maximum reflector tuning efficiency, TuningEfficiencyReflector, at 
lr_min may then be calculated. TuningEfficiencyReflector is typically measured 
in THz/mA. This calculation is performed by finding the first OP that has Ir > 

10 lr_min, for each band except the first band. From this OP (OP[band][col]) and 
the previous OP (OP[band][col-1]), the tuning efficiency may be calculated as 
the change in frequency relative to the change in reflector current. Ir. If there is 
no previous OP, the next OP, (OP[band][col+1]) may be used. If neither the 
previous or the next OP exists, then that band may be ignored. The output 

15 value is the maximum reflector tuning efficiency across all bands. 

The relative variation of output power with coupler current, 
RelPowerVariationCoupler, from lc.nu_min to lc.nu_max may be calculated. 
RelPowerVariationCoupler may be presented in dB. This calculation may be 
performed by finding, for each column, the maximum and minimum power for 

20 the operation points that have a coupler current, Ic, that lies between Ic.nu^min 
and lc.nu_max. Columns that do not reach approximately all the way from 
lc.nu_min to lc.nu_max may be ignored. The average ratio of maximum to 
minimum power across all columns may then be calculated. This average ratio 
may be converted to a dB value, to give RelPowerVariationCoupler. 

25 The relative variation of output power with reflector current, 

RelPowerVariationReflector , from lr_min to lr_max may be calculated. 
RelPowerVariationReflector is typically presented in dB. This calculation may 
be performed by finding, for each band, the maximum and minimum power for 
the operation points that have a reflector current, Ir, between lr_min and 



Attorney Docket: 980.1 131 USUI 



Ir_max. Bands that do not reach approximately all the way from lr_min to 
lr_max may be ignored for this calculation. The average ratio of maximum to 
minimum power across all bands is calculated and converted to a dB value to 
give RelPowerVarlationReflector. 

One other calculation Is to find the first and last OP within each band in 
the area bounded by lr_min <= Ir <= lr_max, lc.nu_min <= Ic <= lc.nu_max. 
Between the first and last OP in each band, the number of OP within the band 
are that are evenly distributed across the band with respect to the column 
indices are counted, to give ColumnCount-2. Those OP that yield the minimum 
and maximum power may be identified and stored in an OP list, OPLIst. 

The next step, step 514, is to perform a phase scan. In this 
measurement, the output power and frequency are measured for each OP in 
OPLIst as a function of phase cun-ent, Ip. A typical result of current and 
frequency measurement is presented in FIG. 19. Curve 1902 shows the 
variation of frequency with Ip, while curve 1904 shows the variation of power 
with Ip. 

One approach to analyzing the phase current data is as follows. First, 
the phase currents associated with phase tuning of 2-n- and 4tt are calculated 
as lp_2TT and Ip_4tt respectively, for all operating points In OPList. For each 
OP, this is done by first finding the start frequency, at Ip = 0, finding the next 
negative frequency hop, typically one that is larger than 0.01 THz, and then 
finding the phase current, lp_2Tr that yields the same frequency as the start 
frequency. This may be found by linear interpolation. Another, similar step 
may be used to find lp_47T. 

The value of lp_max may be calculated as the maximum value, across 
all OP of (lp_2TT + lp_4-rT)/2. For all of the OP in OPList, the actual output 
power for Ip = 0, lp_2Tr and lp_47T may be calculated by scaling the measured 
power values by dividing by FiberCouplingEfficiency. 
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The relative variation of output power with phase current, 
RelPowerVariationPhase , from 0 to Ip_2tt, may be calculated by calculating 
the ratio of the output power for Ip =0 and Ip = lp_2Tr for all OP in OPList, and 
then calculating the average. This average may be converted to dB. 
5 The tuning efficiencies may be calculated for all of the OP in the OPList, 

for the different phase currents, Ip, Ip_2tt and lp_47T. This may be done by 
finding two measurement points that straddle the particular phase current 
value. The tuning efficiency (dnu/dip) is calculated as the ratio of the change in 
frequency with the change in current. 
10 The tuning efficiency, (dnu/dlp)ip_min, at the minimum specified phase 

current, may then be calculated. This may be done using the expression 
(dnu/dlp)ip_min = 1/sqrt(A^ + 2.B.lp_min/Deltanu_GCPRMode), where A = 
1/(dnu/dlp)o, B = 1/(dnu/dlp)ip 2Tr- 1/(dnu/dlp)o, and GCPRMode is the cavity 
mode separation for the entire laser. This value of (dnu/dlp)ip min is 
1 5 approximately equal to the slope of the nu(lp) curve at lp_min. 

Another parameter, TuningEfficiencyPhase, typically measured in 
THz/mA, may be found by taking the maximum value of (dnu/dlp)ip_min. 

The highest output power, P_max, may be found by taking the highest 
power at Ip = 0 and multiplying with 10'^(-0.05*RelPowerVariationPhase). The 
20 multiplication factor of 0.05 originates from the fact that the power is measured 
at Ip = 0mA, whereas an engineering specification may imply that the power 
should be measured at lp_min. It is assumed that lp_min corresponds to a 
maximum tuning of half the cavity mode spacing Deltanu_GCPRMode, in other 
words a phase tuning of tt. Hence, the factor 0.05= 0.5/10. 
25 The lowest output power, P_min, may be found by taking the lowest 

power at Ip = Ip_2tt and multiplying with 10'^(-0.05*RelPowerVariationPhase). 
The overall relative power variation, P_var, may then be calculated as 
P_var = 10*log( P_max / P_min ). 
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New operating points may be created by replacing Ip, which was zero for 
the first OPList, with Ip_2tt, for all OP in OPList. A new OPList may then be 
generated that contains both the old and the new OP. 

The final step 516 in the laser characterization process is the 
5 measurement and analysis of the laser threshold. The measurement may be 
made by measuring the output power, P, as a function of gain current, Ig, for 
each of the operating points generated in step 514. A typical result is illustrated 
in the L-l curve, shown in FIG. 20. 

The L-l curve may then be analyzed to produce a value of the laser 
10 threshold and the differential efficiency. It will be appreciated that different 
approaches may be followed to find these parameters. One such approach is 
now described. 

For each OP in the OPList, the power values are first scaled by dividing 
by FiberCouplingEfficiency, and the maximum output power, P_max, is 

15 determined. For the example presented in FIG. 20, P_max is around 4.6 mW. 
That part of the L-l curve between two thresholds is selected. The first of the 
two thresholds, PJow, may be taken as the maximum value of 0.15 mW and 
0.01 X P_max. The second threshold, P_high is taken as 0.1 x P_max. If there 
are at least a selected number of measurement points on the selected part of 

20 the P(lg)-curve, for example five measurement points, then a straight line may 
be fitted to these points. The threshold current, Ith, Is calculated as the 
intercept of this straight line on the x-axis. The differential efficiency, n, is given 
as the slope of the straight line. 

The minimum and maximum threshold currents, lth_min and 

25 lth_max over all the operation points may then be determined, and the relative 
threshold variation lth_max/lth_min may be calculated. Also, the operating 
points having minimum and maximum differential efficiency, eta_min and 
eta_max, may be identified. 
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The ratio of the highest power of all OP at the minimum specified gain 
current, lg_min, to the lowest power of all OP at the maximum specified gain 
current, lg_max, may be calculated as RelPowerVariation, expressed in dB. 
This is a measure for the maximum output power variation of the laser across 
5 the tuning band, if we allow the gain current to vary between lg_min and 

lg_max. If RelPowerVariation is negative, full equalisation of the output power 
across the tuning band may be possible. 

It will be appreciated that a procedure for characterizing a laser need not 
include all the steps listed herein, or may contain modifications of such steps. 

10 For example, extrapolations and interpolations may be made using techniques 
other than other linear extrapolation and linear interpolation. 

As noted above, the present invention is applicable to characterizing 
laser diodes, and is believed to be particularly useful for characterizing widely 
tunable laser diodes that can be tuned to many different operating modes. 

15 The present invention should not be considered limited to the particular 
examples described above, but rather should be understood to cover all 
aspects of the invention as fairly set out in the attached claims. Various 
modifications, equivalent processes, as well as numerous structures to which 
the present invention may be applicable will be readily apparent to those of skill 

20 in the art to which the present invention is directed upon review of the present 
specification. The claims are intended to cover such modifications and devices. 
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